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Abstract 



We discuss some details of the model proposed in Ref. lV\, in which neutrino oscillation, dark 
matter, and baryon asymmetry of the Universe would be simultaneously explained by the TeV-scale 



D , physics without introducing very high mass scales. An exact discrete Z2 symmetry is introduced, 

under which new particle contents (a real singlet scalar field, a pair of charged singlet scalar 
^ ' fields and TeV-scale right-handed neutrinos) are assigned to have odd quantum number, whereas 

a^ . 

CN| . ordinary gauge fields, quarks and leptons, and two Higgs doublets are even. Tiny neutrino masses 
00 

I are generated at the three loop level due to the exact Z2 symmetry, by which stability of the dark 

■ matter candidate is also guaranteed. The extra Higgs doublet is required not only for the tiny 

. neutrino masses but also for successful electroweak baryogenesis. We discuss phenomenological 

. ^ properties of the model, and find that there are successful scenarios in which above three problems 

X 

^ ' are solved simultaneously under the constraint from current experimental data. We then discuss 

predictions in such scenarios at ongoing and future experiments. It turns out that the model 
provides discriminative predictions especially in Higgs physics and dark matter physics, so that it 
is testable in near future. 

PACS numbers: 14.60.Pq, 14.60. St, 14.80.Cp, 12.60.Fr [August 16, 2009] 



* Address after April 2009: Department of Architecture and Building Engineering, Hokkai-Gakuen Univer- 
sity, Sapporo 062-8605, Japan 



1 



I. INTRODUCTION 



While the standard model (SM) for elementary particles has been successful in describing 
high energy phenomena at colliders, today we have definite motivation to consider a model 
beyond the SM. First of all, observed data for neutrino oscillation indicate that neutrinos 
have tiny masses and mix with each other j^. Second, cosmological data have revealed 
that the density of dark energy and dark matter (DM) in the Universe dominates that 
of baryonic matter [s!]. The essence of DM would be weakly interacting massive particles 
(WIMPs). Finally, asymmetry of matter and anti- matter in our Universe has been addressed 
for a long time as a serious problem regarding existence of ourselves They are all beyond 
the scope of the SM, so that a new model is required to explain these phenomena. 

A simple scenario to generate tiny neutrino masses {m,y) would be based on the see- 
saw mechanism with heavy right-handed (RH) Majorana neutrinos [s^. Then, m,y can be 
described as mj, ~ m'jy/Mfi, where M/j (~ 10^^~^^ GeV) is the Majorana mass of RH 
neutrinos and m^, is the Dirac mass of at most the electroweak scale. This scenario would 
be compatible with the framework with large mass scales like grand unification. The heavy 
RH Majorana neutrinos would generate lepton number asymmetry in their CP violating 
decays that could be transfered into baryon asymmetry [6]. In a supersymmetric model 
with such heavy RH neutrinos, one might find a DM candidate of lightest supersymmetric 
particle H. I„tr„duc«o,> of such large scales, however, causes a problem of h.crarch. In 
addition, the decoupling theorem [8| makes it far from experimental tests. 

An alternative approach is a quantum mechanical generation of small neutrino masses. 
The original idea of radiatively generating neutrino masses due to TeV-scale physics has 
been proposed by Zee ^. Adding an isospin doublet scalar field and a charged singlet 
field to the SM, the neutrino masses are generated at the one-loop level. Phenomenology 
in the Zee model has been studied in Ref. 



10 



Although the original Zee model has 



been excluded by the neutrino data, lots of studies for its extensions and variations have 
been proposed [12] ■ Another possibility for generating neutrino masses via the new scalar 



particles is the Zee-Babu model 13|, , where the neutrino masses arise at the two- loop 



level. Some extensions of these models are discussed for the purpose of the explanation for 
DM and baryon asymmetry [isl . . The extension with a TeV-scale RH neutrino has been 
discussed in Ref. [171], where the neutrino masses are generated at the three-loop due to the 
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exact Z2 symmetry which forbids the Dirac neutrino mass at the tree level unlike the tree- 
level seesaw mechanism mentioned above, and the Z^-odd RH neutrino is a candidate of DM. 
This model has been extended with two RH neutrinos for the description of the neutrino 
data [18|, which however cannot include a mechanism of baryon number generation. The 
idea of simultaneous explanation for radiative generation of neutrino masses and stability of 
DM by introducing various exact discrete symmetries with right-handed neutrinos has been 
used in several models 19|]. Some other models including mechanisms for baryon asymmetry 



20l | and also in the scenario of electroweak 



have been considered in the leptogenesis scenario 
baryogenesis 

In this paper, we investigate some details of a model proposed in Ref. in which 
neutrino oscillation, origin of DM and baryon asymmetry would be simultaneously explained 
by the TeV-scale physics. As we try to built a renormalizable TeV scale model to explain 
these phenomena simultaneously without fine-tuning, we do not impose unnatural hierarchy 
among the mass scales. The model contains an extended Higgs sector with TeV-scale RH 
neutrinos in addition to the SM particle contents. Tiny masses of left-handed (LH) neutrinos 
are generated at the three-loop level under an exact Z2 symmetry. The lightest neutral Z2- 
odd state is a candidate of DM. Baryon asymmetry can be generated at the electroweak 
phase transition (EWPT) by the non-decoupling property [22| and additional CP violating 



phases in the Higgs sector 23|, |2J]. In this framework, a successful model can be built 



without contradiction of the current data. Notice that in this model we do not intend to 
solve the so-called hierarchy problem. The model predicts quadratic divergences in quantum 
corrections to masses of the scalar bosons as in the SM. This model, therefore, has to be 
considered as an effective theory below the cut-off scale of a more fundamental theory around 
at most 10 TeV, below which the self-coupling constants for additional scalar bosons do not 
become larger than the acceptable values for a perturbation calculation. 

We show that there are several possible scenarios, in which the data for neutrinos, lepton 
flavor violation as well as the WMAP data are satisfied. The collider data from experiments 
at LEP, Fermilab Tevatron and B factories at KEK and SLAC are also taken into account. In 
add,tio„ to the o.gi.a. scenano discussed ,„ Ref Q. we .ne.tion another seenano whe.e the 
recent data from PAMELA and ATIC would also be included with the relatively heavy DM 
candidate. Furthermore, we also discuss the other scenario where the anomaly in the data 
from DAMA/LIBRA would be explained by a light DM candidate whose mass is around 5 



3 



GeV. It turns out that in these scenarios all the masses of additional physical particles are 
between (9(10) GeV and 0{1) TeV. We find that the model has discriminative features in 
Higgs phenomenology, lepton flavor physics and DM physics, so that it is testable at current 
and future experiments. 

In Sec. II, we define the model with introducing new particle entries and symmetries, and 
discuss the physical states and their masses and coupling constants. In Sec. Ill, we calculate 
the neutrino mass matrix in our model, which is induced at the three loop level, and discuss 
parameter sets in which all the neutrino data are reproduced under current experimental 
bounds. Sec. IV is devoted to the discussion on the possibility that the lightest Z2 odd scalar 
boson is a candidate of DM. The thermal relic abundance is evaluated in several parameter 
sets, and implication for the physics at direct and indirect search experiments are also dis- 
cussed. In Sec. V, we study the allowed region where the strong first order electroweak phase 
transition is realized. This is required for a successful scenario of electroweak baryogenesis. 
In Sec. VI, we summarize the constraints from the current experimental data on the model, 
and discuss phenomenological predictions at ongoing and future collider experiments and at 
direct/indirect searches for DM. Some formulas are summarized in the Appendices. 



II. MODEL 



A. Particle contents and symmetries 

In addition to the known SM fields, new particle entries in our model are 

^^,^2,S^,r|,N^, (1) 

where $1, $2 are scalar isospin doublet fields with the hypcrcharge 1/2, arc charged 
isospin singlet fields, is a real scalar singlet field, and is the a-th generation isospin- 
singlet RH neutrinos. It turns out that at least two generations are necessary for to 
reproduce the neutrino data. In the following, we mainly consider the minimum model 
with two generation {a = 1,2). We give a comment for the case of more than three 
generations later. We here only note that the mass scale of new particles derived in the 
following sections is not so sensitive against the number of Nji generation. This means that 
the model has high predictability for the mass scale of new particles. 
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$1 


$2 






Z2 (exact) 


+ 


+ 


+ 


+ 


+ 


+ 


+ 






Z2 (softly broken) 


+ 






+ 


+ 


+ 




+ 


- + 



TABLE I: Particle properties under the discrete symmetries. 



In order to generate tiny neutrino masses in the three loop level, in other words to forbid 
the tree-level Dirac neutrino mass term, and at the same time in order to have a stable 
DM candidate, we i,„pose a new (exact) .yn,n,et.y as ,n Ref. fl. wh.eh we .efe. to as 
Z2. We assign the Z2 odd charge to A^^^^, and t], while ordinary gauge fields, quarks and 
leptons and Higgs doublets are Z2 even. 

It has been well known that introduction of the extra Higgs doublet causes a problem 
of dangerous fiavor changing neutral current (FCNC). To avoid FCNC in a natural way, 
we further introduce the (softly-broken) Z2 symmetry in the model, which can be softly 



broken 25j in the Higgs potential. We refer to this symmetry as the Z2 symmetry in 
this paper. Under the Z2 symmetry, there can be four independent types of the Yukawa 
interaction, depending on the assignment of Z2 charges for quarks and leptons. From a 
phenomenological reason discussed later, we employ the so-called Type-X Yukawa coupling 
where charges are assigned such that only $1 couples to leptons whereas $2 does to 



quarks 
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27 



28 



29| . We summarize the particle properties under Z2 and Z2 in TABLE [H 



Under these discrete symmetries, the Yukawa interaction is given by 



J-'Y = -ViiL $i4 - yuiQ^2U^R - VdiQ ^2d^R + h.c 



(2) 



where (L*) is the ordinary i-th generation LH quark (lepton) doublet, and and d^j^ 
(e^) are RH-singlet up- and down- type quarks (charged leptons), respectively. Notice that 



26 



m 



29| is different from that in the 



the Type-X Yukawa coupling defined in Eq. (] 

Type-I or Type-II two-Higgs-doublet model (THDM) [30|. Our Z2 charge assignment for 
quarks is the same as that in Type I, but that for leptons is the same as of Type II. The 
charged Higgs boson couplings to leptons are multiplied by tan/3, while those to quarks are 
by cot j3 in a universal way, where tan/5 = ($2)/('^i)- Therefore, phenomenology of the 
Higgs sector in our model is completely different from that in Type I and Type II. 
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The scalar potential is then given by 



V = - ^^l\^2? - (/i?2*I*2 + h.c.) 

+Ai|$i|^ + A2|$2|' + A3|$in$2r + A4|$I$2r + |y($I$2)' + h.C.j 



+^i\s\' + K\s\' + + Ay + ei^r^ + E {p'^\^'^\ 

a=l ^ 



+ 5^ {«:e,,TOt$,S-r/ + h.c.}, (3) 

a,b=l 

where eab is the anti-symmetric tensor with ei2 = 1- The mass term and the interaction for 
N'^ are given by 

^y = il {\^N-J^'N% - hfie^^N^S^ + h.C.j . (4) 

The parameters /i^25 -^5 ^'^d ^is well as /i" are generally complex. The phases of A5 
and K can be eliminated by re-phasing and $1. The remaining phase of ^"(2 is physical 
and causes CP violation in the Higgs sector, which is necessary for generating baryon 



asymmetry at the EWPT [23|, |2J]. Although the CP violating phase is crucial for successful 
baryogenesis, it does not much affect in the following discussions on neutrino masses, DM 
and the strong first order EWPT required for electroweak baryogenesis. Thus, in the 
following, we neglect the phase of /i^2 (and h"') for simplicity. We later give a comment on 
how the phenomenology could change with non-zero CP- violating phase. 



B. Higgs states and coupling constants 

Because Z2 is exact, Z2 even and odd fields cannot mix. The Higgs doublet fields $j 
{i = 1, 2) can be parameterized as 

VTl^^* + hi + iZi) 

where Vi are vacuum expectation values and satisfy a/w 1 + f | = f 246 GeV), and tan f3 = 
V2/V1. The Z2 even states are diagonalized in mass as a usual THDM by introducing the 
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mixing angles ol and /3, where a is that between CP even neutral Higgs bosons [3 

cos /5 — sin /5 
sin (5 cos (5 

and 







wf 











2;i 




COS /5 — sin /3 
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Z2 




sin P cos /? 




/I 



(6) 







h2 







H 




h 



(7) 



cos a — sm a 
sin a cos a 

where w^ and z are Nambu-Goldstone bosons eaten by and Zj^, and A are charged 
and CP-odd scalar states. CP-odd state can be a mass eigenstate only when CP is conserved 
in the Higgs sector. When CP is conserved h and H are mass eigenstates. Consequently, 
the Z2 even physical scalar states are two CP-even {h and H), a CP-odd (A) and charged 
(if ^) states, like in usual THDMs. Throughout this paper, sin(/3 — a) = 1 is taken, under 
which we define h and H such that h is the SM-like Higgs boson ^. 
The mass formulas for the scalar fields are calculated as^ 



sin2(a - p)M^^ + sin 2 (a - p)M^2 + cos^(a - p)M, 



221 



m% = cos^(a - /5)Mii + sin 2 (a - (3)M^2 + sin^(« - /?)^22, 

2 Ji/r2 A4 -I- A5 2 
m^± = M V , 



mi 



.2 I /„ „„„2 



"^5± = yU5 + (Picos /9-hp2sin 



-h(cricos /3-FcT2sin 



where the mass matrix elements for the CP-even bosons are given by 
= 2(Aicos^/3 + A2sin^/3 + Acos2/3sin2/5)^;^ 



M; 



12 



M 



21 



-2Ai cos^ /3 + 2A2 sin^ /5 + A cos 2/3) sin (3 cos 



M; 



22 



+ -(Ai + A2 - A)(l - cos4/5)t;^ 



(8) 
(9) 
(10) 

(11) 
(12) 

(13) 

(14) 
(15) 
(16) 



where A = A3-I- A4-I-A5. The parameter M relates to the invaiant mass scale fiu in the Higgs 
potential by M = |/ii2|/(sin/5 cos/3) and has the meanings of the soft breaking parameter 
for the Z2 symmetry. 



^ Notice that this is different from the definition that h is always the lighter CP even Higgs boson when the 

SM-like Higgs is heavier than the other one. 
^ For the expressions in the case of sin(/3 — a) = 1, see Appendix A. 
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In terms of mass eigenstates, the scalar interaction terms are given by 



-(pi cos^ (3 + P2 sin^ /3) (^-uj+ + S+S- - (pi sin^ (3 + cos' /3) (^if+if- + 5+5" 

— sin/5cos/5(-pi + p2)(c<;+if- + + 

— (pi cos a cos P + P2 sin a sin P)vHS^S^ — (— pi sin a cos P + P2 cos a sin P)vhS~^S~ 
— -(pi cos' a + p2 sin' a)HHS^S^ — -(pi sin' a + p2 cos' a)hhS^S~ 

— cos a sin a{—pi + p2)hHS^S^ , (17) 



((Ti cos' /9 + (72 sin' P) ^ uj^ + y ^ 



2\ ^2 



2 (di sin^ P + a2 cos' /?) /f+if - 2 J 2 

2 



r /_„:„2^, 2 I rr+rr- . '^^\ 



-sin/?cos/?(-(Ti + (T2)(t^ +a;"if + + zA) — 

77' 77' 
— (o"i COS a cos /5 + cr2 sin a sin P)vH— — (— o"i sin a cos /3 + cr2 cos a sin P)vh— 



1 

■2' 



— -(cTi cos' a + (72 sin' a)HH^ — 2^'^^ ^^^^ + cos' a)hh%^ 



— cos a sin a(— cri + a2)hH—, 



(18) 



and 



£^ = -v/2k [vH^ - sin(a - /5)i/cj+ - cos(a - p)hiu^ 
+ cos(a - P)HH+ - sm{a - P)hH+] S~r] + h.c. . 



(19) 



The Yukawa interactions with quarks and leptons in Eq. ([2]) can then be written as 



^Quarks 



-m,,iU U 



V sin P 



V sin P 



-m^M - ^^^^Hd^d^ - "^^^^hd^d^ 



V sin P 



V sin P 



+1h^zuH-f^u' - ^ cot pAuH-i^u' + "^zdH-i^d' - ^ cot pAdH^i^d' 



V 

+uj~d' 
+H-d' 



V V 

V^m^^ f 1 + 75 \ _ V2"^d' f 1 -75 I j^Y^^ 
V \ 2 J V \ 2 

V^m„i cot P / 1 + 75 \ V^m^t cot P / 1 — 75 



u' + h.c. . (20) 
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/'Leptons 



—m£i 



^Plrrii^ _ . / 1 + 75 



-V 



V 



•i, ,+ 



^J2m(i tan /? _ ■ / 1 + 75 
V 



vfip. cos a 



Bee + 



+ + h.c. 



Z75^' + ^tan/3Af Z75f . (21) 



V cos/? f cos/5 

The coupling constants are constrained by the conditions of theoretical consistencies such 
as perturbative unitarity 



scale A of the model 
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311] or vacuum stability and triviality as a function of the cutoff 
32l. We will take into account these conditions and work in the 



allowed parameter region for A > 10 TeV. 



III. NEUTRINO MASS AND MIXING 



In this section we calculate the mass matrix in our model, which is induced at the three 
loop level, and study the parameter regions which satisfy the current data. The mass term 
for LH neutrinos 

= l^M^^UL^ (22) 
is naturally generated from dimension five operators 

= ^'^vl'H. (23) 

in the low energy effective theory, where represents the neutral component of the Higgs 
doublet, are dimensionless coefficients, and A^v is a dimensionful scale of new dynamics. 
In ordinary models based on the seesaw mechanism, the tiny neutrino mass is essentially 
reahzed by taking a very large scale for A^r jsl: A^v = 10^'^"^^ GeV is required for = C(l) 
to obtain the tiny mass scale comparable to the neutrino data. When A^v is at most TeV 
scales, very small values for are required to describe the data. It would be possible that 
such small would be generated at loop level without making fine tuning on the coupling 
constants in the Lagrangian. 



A. Evaluation of the three-loop induced neutrino mass matrix 

In our model, the LH neutrino mass matrix M^- is generated by the three-loop diagrams 
in FIG. [H The absence of lower order loop contributions is guaranteed by Z2. The charged 
Higgs boson from the doublets and the charged leptons e\ play a crucial role to connect 
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FIG. 1: The diagrams for generating tiny neutrino masses. 

LH neutrinos with the one-loop (box) diagram by the Z2-odd particles. The resulting LH 
neutrino mass matrix is obtained as 

2 

Ml'j = J2CriF{mH±,rns,mN^,m,), (24) 

Q=l 

where 

C^. = AK'tan'(3iy^yt)iyrh';), (25) 



and the loop integral function F is given by 



oo 



X xdx{Bi{-x,mH±,ms) - Bi{-x,0,ms)} — ■ ^ ■ -], (26) 



2 



with TTif representing the mass of the field / and yf^ = \/2meJv. The function B 



is the 



tensor coefficient in the formalism by Passarino-Veltman for one-loop integrals 33|. The 
detailed calculation of F is shown in Appendix B. 

In Fig. [21 we show the magnitude of the integral function F as a function of m^y for 
several values of m^± and m^. It can be seen that F becomes smaller for larger values of 
m^± . For < 500 GeV F decreases monotonically as rrijsf^ grows, and for greater values 
of it mildly increases but finally turns to decrease as itLj^^ grows. The dependences on 
are however not very sensitive. Numerically, the magnitude of F is of order lO'* eV 
in the wide range of parameter regions of our interest. We note that the dependence on 
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F(xlO eV) 




FIG. 2: The contour plot of the values F as a function of and for = 50 GeV and 
m^± = 100 GeV. 

the mass of the charged Higgs boson is quahtatively similar to that on m^±, so that F 
becomes smaller for larger values of m^^. 

Magnitudes of k and tan (3 as well as Fijn^^mg^mj^^^mj^) determine the universal 
scale of M"j, whereas variation of /i" {i = e, /i, r; a =1-2) reproduces the mixing pattern 
indicated by the neutrino data j^. 



B. Parameters to reproduce the data of neutrino masses and mixings 

The generated mass matrix M^j in Eq. fl24|) of LH neutrinos can be related to the neutrino 
oscillation data by 

M^. = f/,,(M,^i,g),,(f/^),„ (27) 

where M^^^^g = diag(mi, m2, ma). For the case of the normal hierarchy we identify the mass 
eigenvalues as mi = 0, m2 = \/^"^iar ^^'^ "^3 ~ V ^"^atm; while for inverted hierarchy 
mi = a/ Am^tjn, m2 = a/ Am^^m + ^"^soiar ^^"^ ms = are taken. The Maki-Nakagawa- 
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Sakata matrix 



34| is parameterized as 



U 



1 










'^13 





5i3e^' 




•^12 


■^12 







1 











^23 


"^23 







1 







^■^12 


Cl2 


















^■^23 


•^23 






' 


Cl3 










1 













(28) 



where Sij and Cij represent sin 6ij and cos % respectively with 6ij to be the neutrino mixing 
angle between the ith and jth generations, and 6 is the Dirac phase while a and f3 are 
Majorana phases. For simplicity, we neglect the effects of these phases in the following 
analysis. Current neutrino oscillation data give the following values 



Am 



solar 



7.65 X 10-^ eV^ 



l^"^atml - 



2.4 X 10-^ 



sin^ 9i2 ^ 0.3 , sin^ ^23 ^ 0.5 , 



sin^^ia < 0.04 



(29) 
(30) 



In order to study parameter sets that satisfy the current experimental data, we need to 
discuss constraints from the /i — > 67 results. This process is induced by one loop diagram 
of a right handed neutrino Nji and a charged scalar boson through Yukawa coupling 
[i = e and fi). The branching ratio is given by 



5(/i 67) 



where Fo(x) 




(31) 



Qx + 3x + 2x — Qx — \nx)/Q(l — x) . Comparing this to the current 



experimental bound, B{fi 67) < 1.2 x 10^^^ 35], the mass of Nj^ or is strongly 
constrained from below. 

Under the natwra/ requirement that h'^ = 0{1), and taking into account the fi ^ ej search 
results y, we find that m^. > 0(1) TeV, m^± = 100-200 GeV and Ktan/5 = C(10). In 
addition, with the LEP direct search data for Higgs bosons and precision measurement 
data [2], possible values for the scalar masses uniquely turn out to be fn^j^ ~ 171^^ ~ 100 



GeV and ~ (9(100) GeV for sin(/5— a) = 1. As we see in later sections, relatively heavier 
S"^ (^5± ^ 400 GeV) is favored from the discussion on DM and electroweak baryogenesis. 
It is known that in the Type X Yukawa interaction a light ('^//i ^ 300 GeV) is not 
excluded by the b —>■ s'j data [36^. This is the reason that the Yukawa coupling in Eq. ([2]) 
has been employed in our model. 

For values of ((i = fi,T)), we only require that hfyi ~ C(?/e) ~ 10~^, since we cannot 
avoid to include the hierarchy among yf^. Several sets for /i" which satisfy the neutrino 
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Set 


Mass (TeV) 




Coupling constants 


LFV 


(hierarchy, sin^ 2^x3 ) 


niri 


nis 




Ktan (3 












^2 




A (normal, 0) 


0.05 


0.4 


3 


29 


2.0 


2.0 


0.041 


-0.020 


0.0012 


-0.0025 


6.8 X 10"^^ 


B (normal, U.14j 


0.05 


0.4 


3 


34 


2.2 


2.1 


0.0087 


0.037 


-0.0010 


0.0021 


5.3 X 10^^^ 


C (inverted, 0) 


0.05 


0.4 


3 


66 


3.8 


3.7 


0.013 


-0.013 


-0.00080 


0.00080 


4.2x10^12 


D (inverted, 0.14) 


0.05 


0.4 


3 


66 


3.7 


3.7 


-0.016 


0.011 


0.00064 


-0.00096 


4.2x10^^2 



TABLE II: Values of hf which satisfy neutrino data and the constraint from /i ^ 67 for m,j = 50 



GeV, Tn^^- 



lOOGeV and m^yi = 1x1^2 . 



Set 


Mass (TeV) 




Yukawa couplings 


LFV 


(hierarchy, sin^ 2^13) 




ms 


" R 


Ktan/3 


hi 


hi 


hi 


K 


hi 


hi 




E (normal, 0) 


0.7 


0.9 


3 


33 


3.0 


2.0 


-0.014 


0.057 


-0.0028 


0.0021 


7.2 xl0~^2 


F (normal, 0) 


0.005 


0.4 


3 


29 


2.0 


2.0 


-0.019 


0.041 


-0.0024 


0.0012 


6.6 xl0"i2 



TABLE III: Values of hf which satisfy neutrino data and the constraint from /i ^ 67 for 



700 



GeV (Set E) and = 5 GeV (Set F). The other parameters are taken to be rn^j^ = lOOGeV and 
m^2 for the normal hierarchy. 



data are shown in TABLE [TTl with the predictions on the branching ratio of /i ^ 67 for 
rriri = 50 GeV, mH± = 100 GeV and ttlj^i^ = = 3 TeV. Set A and Set B are rather 
standard choices in our model assuming the normal hierarchy, while C and D are those 
assuming the inverted hierarchy. We note that larger values are required for n tan (3 in the 
case of the inverted hierarchy. On the other hand, small tan (3 is favored by the realization 
of electroweak baryogenesis as we will discuss later^. Therefore, our model turns out to be 
better compatible with the normal hierarchy scenario [37]. 

Set E and Set F are chosen to potentially include the results from PAMELA (and ATIC) 
and those from DAMA, respectively, as discussed later. The parameters describing the data 
are listed in Table IIIII 



In the THDM, tan (3 > 10 would not be favored for successful electroweak baryogenesis under the EDM 
constraints 



24|. Although we do not discuss CP violating effects in our present analysis, we keep in mind 



this constraint and just impose the condition that tan/? is not larger than about 10. 
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We also examined the case with three generations for RH neutrinos. Unhke the case with 
two generations, the mass of the smallest mass eigenvalue for LH neutrinos can be nonzero. 
It is however found that larger values for /ttan/? are prefered to reproduce the neutrino 
data, so that it would be rather unnatural as compared to the case with two generations. 
Therefore, we concentrate on the model with two generation RH neutrinos in the rest of the 
paper. 

In Fig. [3l the contour plots for the branching ratio of /i ^ 67 as a function of fn^q^ and 
for values of /i" to be of Set A - Set F in Table II and III at = 100 GeV. The 

dependence of the branching ratio on the mass of 5*^ is rather mild. On the other hand, the 
branching ratio is strongly dependent on the mass of N^, so that the lower bound on ttL]^^ 
is obtained for a mass of 5^. 

IV. DARK MATTER 

Since Z2 is exact, the lightest Z2-odd particle is stable and can be a candidate of DM 
if it is neutral. In our model, N'^ must be heavy (m^^ ~ a few TeV), so that the DM 
candidate is identified as 77. Since 77 is a singlet under the SM gauge group, the interactions 
with Z2 even particles are only through the Higgs coupling. It is, however, worth noticing 
the presence of interactions through k, coupling, which is absent in the simplest gauge singlet 
scalar DM U. 




A. Thermal relic abundance 

In principle, 77 has various annihilation processes into 77, //, W^W~ and so on. The 
annihilation into a W boson pair or a Higgs pair leads to too rapid annihilation. Thus, 
relatively light 77 would be favored in order to kinematically close these annihilation modes. 
When 7; is lighter than the W boson, r] predominantly annihilates into two photons through 
one-loop diagrams by and S*^, as well as into hh and r+r^ through s-channel Higgs {h 
and H) exchange diagrams: see Fig. HI 
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FIG. 3: The contour plot of ^ 67 branching ratio in the m^y — ms plane for the values of h° 
given in Set A to Set F. 

The relevant annihilation rate is evaluated for sm{(5 — a) ~ 1 as 



1 



IGtt V IQtt 



+ 



IGtt 



(7i cos^ /5 + (72 sin^ P (o"! + cr2) cos^ P 



s - ml + inihTh s - mjj- + irriHTH 



2 1 
smzv 



IGtt 



0"! cos^ /5 + o"2 sin^ P ((Ti — (72) sin^ /3 



I ^ " (32) 

s -ml + iiTLhTh s - mjj + imuTH 

where and Th are the decay widths of h and H respectively, and s is a usual Mandelstam 
variable. The second and third term correspond to the annihilation into fermion and anti- 
fermion. The annihilation into two gammas is given by the first term whose expression was 
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FIG. 4: Feynman diagrams for the DM annihilation. For the one- loop diagrams for 7777 — > 77 only 
those which contain the coupling constant k are shown because they give dominant contributions 
in our parameter choices. 

obtained under the s-wave approximation. From their thermal averaged annihilation rate 
(crv), the relic mass density fi^/i^ is evaluated for 7777 — > 77, bb and r+r~ as 



= 1.1 X 10^ 



GeV-\ (33) 

7d 



where Mp is the Planck scale, g^, is the total number of relativistic degrees of freedom in the 



thermal bath, and is the decoupling temperature [391]. 

Except for the region near the resonance annihilation by s-channel Higgs (h and H) 
exchange where ~ mh/2 and mn/^, the DM relic density is basically determined by the 
annihilation into 77 whose cross section essentially depends on only two free parameters; k 
and m^j^. Since we already know n = 0{1) from neutrino masses discussed in the previous 
section, we obtain m^^ > 400 GeV from the r] DM abundance result. 

FIG. [5] shows Qrjh"^ as the function of for sin(/5 — a) = 1. The parameters are 
chosen as (ai, (T2, k tan/?, tti/^, mj^) = (0.05, 0.03, 30, 120GeV, lOOGeV). Strong annihilation 
can be seen near 50 GeV ~ ^h/'^ (60 GeV ~ wi/i/2) due to the resonance of H (h) 
mediation. The VLr^h?_ is proportional to for the dominant annihilation into 77. The 
data (l^DM^^ ~ 0.11 [3]) indicate that is around 50 — 65 GeV for m^± of 400 GeV. If S 
is much lighter than 400 GeV, the resultant Vtrjh? is below 0.1 for most of range. For a 
heavier m^±, only the resonance annihilation would provide the desired VLrfh?. 
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FIG. 5: [Left] The thermal rehc abundance of ry and contributions from tree and one-loop processes. 
[Right] The thermal relic abundance of r] and its ms± dependence. In both figures, parameters are 
taken as (k, cji, cj2, m^, m^, tan/?) = (3, 0.05, 0.03, 120GeV, lOOGeV, 10). 



B. Constraints from DM direct searches 



40] 



Generally speaking, WIMP DM is detectable by direct DM searches such as CDMS 
Scatterings off of WIMP DM with nuclei take place in two different ways, depending on 
the nature of WIMP Q]. One is spin-dependent scattering whose cross section depends 
on the total spin of target nuclei, while the other is spin-independent (SI), r/ is a scalar 
particle, therefore the relevant process is a SI scattering through t-channel Higgs {h and H) 
\. The SI cross section for a proton is given as [3] 



exchange 



with 



a. 



SI 



mt 



7r{m^ + Trip) 



-f 



A. 



q=u,a,s ^ 



(p) fi I A \^ f ( 

g=c.b,t 



(34) 



(35) 



where rrip and are the proton and each flavor quark masses. The hadronic matrix elements 



f^^ and f^Gq defined as 

T^vfrl = {P\mqqq\p), 
and we adopted the values in Ref. 



c(p) 



(p) 

TGq 



1 



(p) 

Tq^ 



(36) 



411 ] for the following estimation. In our model, the 
coupling constant of the effective interaction between r] and quark, Cint D fqqqvVy is given 
as 



rrir, 



'—(Ti sin a cos /3 + (72 cos a sin /3) cos a (cti cos a cos /5 + (72 sin a sin /3) sin a . 



2ml 



sin P 
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Until now, only null results are reported from all experiments of direct DM search. The 



current most stringent bounds are given by XENON 10 42| and CDMS II as a^^ ^ 5 x 10~^ 
pb. We find the coupling constants (Tj < (9(10^^) are consistent. 

Remember the most of parameter region of the minimal singlet scalar DM is excluded 
by null results of DM direct search and only resonance annihilation regioii, which would 
mean certain fine-tuning between DM mass and Higgs mass, is consistent [38!]. Hence, in 
contrast with this, it is worth emphasizing that such a tuning to realize large resonance is not 
necessarily required in our rj DM, because the one-loop processes give sizable contributions 
to annihilation for ms — 400 GeV. Another singlet scalar DM model without necessity of a 



resonance due to different mechanism is also recently proposed 43 1. 



C. Variations with a particular choice of parameters 

So far, we have discussed properties of rj DM on a standard basis according to parameter 
sets in Table HIl However, there are a few recent observational results which would imply 
that the nature of DM could differ from the standard WIMP. 

The first example is anomalous excesses in the observed positron flux. The PAMELA data 
shows that observed positron flux exceeds the normally expected astrophysical background 
flux for energy range above a few hundred GeV with a rising spectrum in /{^^ + 
) [44]. In addition, ATIC also reported a bump around E ~ 600 GeV in the spectrum |45|. 



46j or decay |47] of DM 



These excesses could be interpreted as a signal from annihilation 
particles. As mentioned above, we have argued that the mass of rj is most likely around 
50 GeV in order to avoid its annihilation into H^-boson or Higgs boson. However, if we 
accept negligibly tiny cTj coupling constants to suppress these annihilation via the s-channel 
Higgs exchange, it is possible to realize that r] is as heavy as several hundreds GeV in Set 
E in Table IIIII Indeed, in such a case, the annihilation into a pair of charged Higgs bosons, 
which subsequently decays dominantly into a tau lepton and a neutrino due to the Type-X 
Yukawa interaction, becomes the dominant mode and the resultant fi^/i^ ~ 0.1 corresponds 
to rj mass of several hundreds GeV The tau lepton produces electrons and positrons as 



well as hadrons. Hence, at least as far as the energy scale is concerned, the ATIC anomaly 



We estimated this from Eq. (|32p . which might be too simple for a heavy 77. However, we suppose that 
the error would be just a factor difference. 
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FIG. 6: A large enough scattering cross section of light DM with a proton to account for the 
DAMA/LIBRA claim. 



might be account for a heavy rj. The detailed calculation, paying attentions to the spectrum, 
will be presented elsewhere j^sl^. 

Next, the results at the DAMA/LIBRA experiment claim a significant annual modulation 
in their DM detection rate which also could be caused by WIMPs [50], whereas all other 
experiments do not find any signal yet. These two results appear to conflict each other, but 
a few compatible scenarios are possible, e.g., inelastic scattering DM |5l|, nevertheless. One 
of them is a light WIMP with the mass of 3 GeV < < 8 GeV and a large scattering cross 
section with nuclei of O{10^^) pb [52] (as in Set F in Table UTTl) . Such a large can be 
available for ^2 ~ 0.4 as seen in FIG. [61 With a slightly heavier charged scalar S to suppress 
one-loop contributions, the desired thermal abundance fi^/i^ ~ 0.1 can be obtained almost 
independently from in the relevant mass range. 



Within the framework of three loop generation of neutrino masses, the explanation of positron excess was 



already examined in Ref. 49| by extending the other model in Refs. 171 Il8l|. 
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V. STRONG FIRST ORDER PHASE TRANSITION FOR SUCCESSFUL ELEC- 
TROWEAK BARYOGENESIS 

The model satisfies the Sakharov's conditions for baryogenesis j4|: baryon number vio- 
lation occurs via sphaleron effect at high temperatures, C and CP violation is automatic 
in the electroweak theory and additional CP violating phases are in the Higgs sector and 
in the Yukawa interaction, and the condition of departure from thermal equilibrium can be 
realized by the strong first order EWPT in the electroweak baryogenesis scenario. In this 
section, we examine the region in the parameter space where the strong first order EWPT 
is realized. 

In this paper, we take the SM-like limit for the mixing between h and H (sin(a— /?) = —1), 
where h is the SM-like Higgs boson while H does not have VEV. In such a case, the one-loop 
(both zero and finite temperature) effective potential is described in terms of a unique order 
parameter (= {h)Y. All extra scalar bosons (if. A, H^, S^, rf) approximately behave 
just like additional particles running in the loop in the effective potential. The following 
relations and equations do not depend on the value of tan/?, so that we apply our result 
here to the case of tan (3 = C(l) - C(10). 



A. Effective potential 

We first consider the one- loop effective potential at zero temperature Ve/zfv^] = Vtreefv^] + 
Ay[(y9]. The one-loop contribution Al^[y9] is given by 

-^^M -^.EV'.f-^fm'f])' {m™^ - 1} , (38) 

where ip = {(j)) = v + {h) , is the color factor, Sf (A^s^) is the spin (degree of freedom) of 
the field / in the loop, M/[y9] is the field dependent mass of /, and Q is an arbitrary scale. 

^ The effect of the two stage phase transition in the case with a large tan /? ^ 1 is neglected in our 
analysis |53j. 
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For sin(/3 — a) — 1, the field dependent mass functions are given by 
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The renormahzed mass of the SM hke Higgs boson h and the one-loop corrected hhh coupling 
are given by the conditions 

d 



d(p 
dip' 

g3 



ip=v 



(p=v 



0, 
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eff 
hhhi 



(45) 
(46) 
(47) 



when sin(/3 — a) = 1. 



B. Electroweak phase transition 



The additional part of the one-loop effective potential at finite temperature is expressed 



by 



where //(a/) is 



poo 

Is{af) ^ dxxMog[l-e-V^], (boson), 
Jo 

poo 

Ipiaf) = / fix x^log[l + e"V^^^], (fermion), 
Jo 



(48) 



(49) 
(50) 
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for the loop particle / with a/ = m{ip)f/T. For a given mass parameter set, the critical 
temperature T^. and the critical expectation value ipc (7^ 0) are obtained as solutions of 



Veff[(fic,Tc] = 0, 

d 

— K//[<^c,T,] = 0. 



(51) 
(52) 



The expression by using the high temperature expansion [5J] (see Appendix D) is useful 
to see the structure in an analytic way. When sin(/3 — a) = 1, m|^,m^,m^± ^ M^, 



m 



|± ^ ;u|± and ^ /i^, the effective potential at finite temperatures Veff[ip,T] 



Vireel'^] + ^^[v^] + '^VtIv^T] cau bc approximately expressed as 



(53) 
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(54) 
(55) 
(56) 



(57) 



where loga^ = 21og47r — 27^ and loga^ = 21og7r — 27^, and 

B - -^^^ {e^t- + 3m| - 12m^ + ^m^ + m\j + m\ + 2m^± + 2m|± + mjj . (58) 
We analytically obtain the critical temperature and the order parameter ip^ at as 



% = To 
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2ET, 



(59) 



In order to satisfy the sphaleron decouphng condition in the broken phase, it is required 



that 



55] 



> 1. 



(60) 
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FIG. 7: The region of strong first order EWPT is shown in the m^-m^-t plane. Shaded (yehow 
colored) area is excluded by the condition in Eq. (|60p . The invariant mass parameters M and fj,^ 
are taken as M = 100 GeV and fig = 200 GeV. 

In terms of high temperature description of the effective potential, Eq. (150]) implies the 
condition 

1^ > 1. (61) 



A large value of E is important for the strongly first-order EWPT [22], for which only the 
bosonic loop effects can contribute. In the case of the SM, the value of the coefficient E 
mainly comes from loop contributions of weak gauge bosons, so that the above relation gives 
the upper bound on which is much below the lower bound from the LEP experiment. 
However, in some new physics models, the situation can be improved due to the contributions 
of additional bosonic fields which couple to h. In our model, there are many additional 
scalars running in the loop so that a larger E can be easily realized. Consequently, the first 
order EWPT is possible without contradiction with the LEP data in a wide region of the 
parameter space. 

In FIG. d the allowed region under the condition of Eq. ( l6Ti) is shown for rrih ~ 120 
GeV, rrifj ~ m^±{p=^ M) ^ 100 GeV and sin(/5 — a) ~ 1. The invariant mass parameters 
M and /i^ in Eqs. ([74]) and ([TSj) are taken as M = 100 GeV and /i^ = 200 GeV. The 
condition is satisfied when m^± > 360 GeV for ~ 100 GeV and when > 340 GeV 

23 



for — 200 GeV. If we take larger values for M (< m^) and jJ-s < "i^i, then the area of the 
allowed regions becomes smaller: in order to obtain similar magnitude of non-decoupling 
effects a heavier A or is favored. The result is not sensitive to tan (3. We note that 
the approximation with high temperatures becomes worse for larger M (M ^ T) but the 
conclusion here is not qualitatively changed for our parameter set even when we evaluate it 
without high temperature expansion. 

VI. PHENOMENOLOGY 

A. Constraints on the parameters from the current data 

Before going to the discussion on the prediction of the model, we here summarize the 
constraints on the model from the following current experimental data and the theoretical 
requirements. 

The // — > 67 results: 

The constraint from the data for — > 67 and the natural generation of the tiny neutrino 
masses require that rrij^cx > 0{1) TeV and ms± is several hundred GeV. 

Tiny neutrino masses with 0{1) coupling constants: 

We require no unnatural fine tuning on the coupling constants k and /i"; i.e. 
K ~ /i^ ~ 0{l). In this case, the data prefer m^± < 100 GeV, mg± = O(IOO) GeV and 
Ktan/3 ~ 30. We do not explain the mass hierarchy among charged leptons, so that we 
consider scales of /if to satisfy 0{h1y}) ~ 10"^ ~ 0{ye). 

The data from the LEP and from the Tevatron: 
Throughout this paper, we take 

m^± = rriff, and sin(/3 — a) = 1. (62) 

In this case, the Higgs potential has the global custodial SU{2)v symmetry. The bound 
from the LEP precision measurement especially those for clcctroweak rho parameter is 
easily satisfied. The scalar h plays the same role as the SM Higgs boson. The other scalar 



24 



bosons do not contribute to the rho parameter. Therefore, we have similar mass bound on 
h to that on the SM Higgs boson from the LEP data; i.e., rrih = 114-160 GeV. The LEP 
lower bound for charged Higgs boson mass should be mH± > 100 GeV, which must also 
be respected. Combining this with the requirement of natural generation of tiny neutrino 
mass, we obtain = rrij^ ~ 100 GeV. A little bit larger values are also allowed when we 
take small fine tuning for the coupling constants. 



The 6 — >• S7 results: 

As discussed in previous section, the results from b ^ S'-f strongly constrain the mass 
of charged Higgs bosons in the Type-II Yukawa interaction in THDMs. It gives the lower 
bound > 295 GeV [56|. In order to avoid this constraint and in order to have a 

compatible value of fri^j^ with the neutrino data, we choose Type-X Yukawa interaction in 
our model. 



The g — 2 results: 

The predicted muon g — 2 in this model with the above parameter sets of new particle 
mass and Yukawa couplings to reproduce neutrino mass is far below the current experimen- 



tal bound 



571]. In addition, the electron g — 2 is also below the current bounds even for 



order unity couplings because it is suppressed by the electron mass [58|. Thus, there are 
no significant constraints on the model from g — 2. 



The dark matter data: 

The DM thermal relic abundance fi^/i^ ~ 0.1 (WMAP) requires < 65 GeV for 
ms > 400 GeV, or ~ ^h/2 or 171^/2. The null result of direct DM searches corresponds 
to couplings a < 10"^ for m,, > 20 GeV. 

Strong first order phase transition: 

We need non-decoupling property in the Higgs sector to realize the strong first order 
phase transition of electroweak symmetry. As, the mass of H and are constrained to 
be around 100 GeV, these fields cannot give such strong non-decoupling loop effects. We 
consider the case where A and are relatively heavy so that the sufficient non-decoupling 
loop effects can be generated. Concretely, we may consider uia = 100 — 200 GeV, and 
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FIG. 8: Branching ratios of the SM hke Higgs boson h for tan/3 = 3 (left) and tan /3 = 10 (right) 
for sin(/? — a) = 1 and = 48 GeV. The couphng constants are taken to be a\ = 0.05, o"2 = 0.03, 
pi = 1.0 and p2 = 1.0. 



m^± > 400 GeV. 



B. Predictions 



We now show phenomenological predictions in the scenario in ( !63|l in order. Here we 
mainly work on the scenario which can simultaneously solve the above three issues under 



the data 



36|: 



sin(/? — a) = 1, 
TUh = 120GeV, 



Ktan/S ~ 30, 

m^j ^ mH±{^ M) ~ lOOGeV, 



niA ~ 100 - 200GeV, m^± ~ 400GeV, 



(63) 



= 40 - 65GeV, 



R 



m^2 = 3TeV. 



By this scenario, the model can explain the origin of tiny neutrino mass and mixing, the 
origin of DM, and the strongly first-order phase transition which is necessary for successful 
electroweak baryogenesis. This can be realized without assuming unnatural hierarchy 
among the coupling constants. All the masses are between (9(100) GeV and 0{1) TeV. 
They are indicated by the data, so that the model has a predictive power. 



Invisible decays of the Higgs boson to a DM pair: 

The scalar boson h is the SM-like Higgs boson, whose coupling constants to the SM 
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FIG. 9: The cross section for the scattering process with a proton at the direct search experiment. 
The upper (lower) curve corresponds to the parameter cj2 = 0.03 (0.01). The other parameters are 
taken to be ruh = 120 GeV, thh = 100 GeV, cJi = 0.05 and tan/? = 10. 

particles coincide with that of the SM Higgs boson at the tree leveL Therefore, its production 
mechanisms are the same as those in the SM. However, when < m/j/2 it decays into a 
DM pair 7777, namely the invisible decay of Higgs as a common prediction of singlet scalar 
DM (such as Higgs portal model) 59| , whose decay rate is given by 



V 



4mH 



327rm^ 



rlo"! cos^ /5+cr2 sin^ P\' 



The branching ratio in our model is shown in Fig. [HI When rrih = 120 GeV and 



(64) 



48 



GeV, the branching ratio oi h ^ rji] amounts to about 36 % (34%) for tan/? = 3 (10). For 
rrih = 120 GeV and = 55 GeV, it is about 25 % (22 %) for tan/3 = 3 (10). For heavier 
h than 120 GeV, the branching ratio oi h ^ rji] is smaller because the mode h WW* 
becomes larger. The decay rate is strongly related to the DM abundance: see Eq. ( I32l) we 
may be able to reconstruct at least the annihilation into //, so that our DM scenario can 
be testable at the Large Hadron Collider (LHC). 



DM searches: 

Our DM candidate 77 is potentially detectable by direct and indirect DM searches. For 
the coupling constants a, ~ 10^^ and the mass ~ 50 GeV, the typical SI cross section 
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is of order of 10 ^ pb and within the reach of near future direct search experiments such 

60l |. In Fig. [9l the SI cross section of DM for a proton a^^ 



as superCDMS and XMASS 



is plotted. In addition, indirect searches by observing cosmic rays which would be a signal 
of DM annihilation is also promising. It is regarded that a line gamma is the smoking gun 
of DM annihilation, rj has a large cross section of annihilation into two photon. Although 
this is a one-loop process, the cross section is not so much suppressed because of a large 
coupling K and light particles inside loops {S^ and H^). As a similar case, one may recall the 
inert Higgs doublet model where the iy+-boson loop gives a sizable contribution 6l(]. The 
FERMI satellite 62], used to be called GLAST, sensitivity with the solid angle of detector 
AQ ~ 10~^ is expected to be > 10~^GeVcm"^s^^ for (9(10) GeV energy. Line gamma- ray 
emissions could be observed for the Navarro, Frenk, and White (NFW) density profile 63] 
with the astrophysics dependent dimensionless function J ~ 10^, because the flux originated 
from rj annihilation is expected to be as large as 



g^^ . 5 X 10- (^) (^^^) J An IGeVc.- V, (65) 

with E = m„. 



The non-decoupling effect of on the hhh coupling: 

As we have seen the condition for the first order electroweak phase transition requires 
non-decoupling property for or A. It is known that such non-decoupling property affects 
physics of the renormalized hhh coupling. In terms of the renormalized mass parameter rrih, 
the one-loop corrected tri-linear coupling is obtained from the effective potential 6J], 




The deviation from the SM prediction 



is defined by AA/j/i/i/A^{{(SM), where AXhhh = — A^{{(SM). The quantum effect on 
the hhh coupling amounts to more than 10-20 % (see FIG. [10]), due to the non-decoupling 
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FIG. 10: The contour plots represent deviations in the hhh coupling from the SM value. For the 
invariant mass parameters, M = 100 GeV and /i^ = 200 GeV are taken. The region of the strongly 
first order phase transition at this parameter choice is also shown. 



property of A and S*^, which is required for successful baryogenesis [22, l2J]. Thus, it should 
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66 



67|. 



be testable at the International Linear Collider (ILC) 

We note that in addition to the non-decoupling effect on the hhh coupling the non- 
decoupling property of the charged singlet scalars S"^ can be further tested by measuring 



77/1 vertices 



Phenomenology of the Type-X Yukawa interaction: 
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291], the phenomenology 



Because of the Type-X Yukawa coupling defined in Eq. ([2]) 
for the extra Higgs bosons is different from that in Type I or Type II especially for tan f3 > 2. 
In our model H and A {H^) can predominantly decay into t~^t~ {t^u), while in the Type 
II THDM the main decay modes of H and A (H^) are bb (rz/) when the masses are around 
100 GeV and sin(/5 - a) = 1. 

Recently, the phenomenology in Type-X THDM has been studied in the similar 
parameter choice in Ref. 28|]. The physics of Type-X Yukawa interaction can be tested at 
the LHC with the low luminosity (30 fb~^) via the single direct (associated) production 
processes gg ^ ^ t~^t~ and fi~^fi~ {pp —>■ bb^ bbT^T~ and bbfi~^fi~), where $ = A and 
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FIG. 11: Production cross section of a S~^S pair via the Drell-Yan process at the LHC {^/s = 14 
TeV) as a function of the mass of S^. 



H . At high luminosity (300 fh'^) , pp ^ W* AH"^ and pp ^ W* ^ HH^ [68 1 can also 
be used to discriminate between the types of Yukawa interaction when the masses are not 
too large. In the minimal supersymmetric standard model (Type II Yukawa interaction) 
the main signal would be the hhrv final state while that in the Type-X the final main states 
would be T^T"T^v or fi^fi^r^u. 



The phenomenology of the charged singlet scalar field S^: 

The physics of Z2-odd charged singlet is important to distinguish this model from 
the other models. At the LHC, they are produced in pair via the Drell-Yan process [h|. In 
Fig- [m the cross section for pp S^S^ at the LHC is shown as a function of the mass 
m^±. The cross section amounts to 0.5 fb for m^±, so that more than a hundred of the S^S~ 
events are produced for the integrated luminosity 300 fb~^. The produced bosons decay 
as S*^ — >■ H'^rj, and if^ mainly decay into t'^u due to the Type-X Yukawa coupling when 



tan/3 > 2 [28|. The signal would be a high-energy hadron pair [69|| with a large missing 
transverse momentum. 

The charged singlet scalar bosons in our model can also be better studied at the 
ILC via e+e~ — » S^S~ shown in Fig. [TWLeft). The total cross sections are shown as a 
function of for several values of the center-of-mass energy y/s in Fig. [T^Left). The 
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FIG. 12: Feynman diagrams for the processes of e+e S^S (Left) and e e S S (Right). 




100 200 300 400 500 -1 -0.5 0.5 1 



rris (GeV) cosG 

FIG. 13: [Left] Production cross sections for e^e~ — > S^S^ via the s-channel gauge boson (7 
and Z) mediation (dotted curve), the t-channel RH-neutrino {N'^) mediation (solid curve), and 
both contributions (dashed curve) for ^/s = 300, 500 and 1000 GeV. The masses for RH neutrinos 
are taken to be = = 3 TeV. [Right] The angular distribution of the e~^e~ — > S~^S~ at 
v/s = 1 TeV for m.,i (= m„2 ) = 0.6, 1, 2 and 3 TeV. In both figures, h], = = 2.0 are taken. 

other relevant parameters are taken as rrij^i = m^2 = 3 TeV and hi = hi = 2.0. Both the 
contributions from the s-channel gauge boson (7 and Z) mediation and the t-channel RH 
neutrino mediation are included in the calculation. The total cross section can amount to 
about 100 fb for = 400 GeV at ^/s = 1 TeV due to the contributions of the t-channel 
RH neutrino- mediation diagrams with 0{1) coupling constants /i" . The signal would be a 
number of energetic tau lepton pairs with large missing energies. Although several processes 
such as e~^e~ W^W~ and e'^e^ H^H^ can give backgrounds for this final state, we 
expect that the signal events can be separated by kinematic cuts. In Fig. [T3l(Right). the 
dependences on the scattering angle in the differential cross section da / d{cos 6) are shown. 
For m,n = m,j2 =3 TeV, the special behavior in the angular distribusion is more insensitive 

R "r 

than the cases with lighter values for m^a. 



31 



100 




1 I 1 , ^ , ^ , ^ 1 

100 200 300 400 500 

ms (GeV) 



FIG. 14: Production cross sections for e^e^ via the t-channel RH- neutrino {N^) media- 

tion (solid curve) for y/s = 300, 500 and 1000 GeV. The masses for RH neutrinos and the relevant 
coupling constants are taken to be m,n = m,,2 =3 TeV and hi = = 2.0, respectively. 

Finally, there is a further advantage in testing our model at the e~e~ collision option of 
the ILC, where the dimension five operator e^e^ S^S^, which appears in the sub-diagram 
of the three loop induced masses of neutrinos in our model, can be directly measured. The 
production cross section for e~e~ S^S^ [t-channel mediation: see Fig. [T2]( Right)] is 
given by 

2 



ax 

aie^e^ — > S^) = I dt- 



(68) 



Because of the structure of our model that the tiny neutrino masses are generated at the 
three loop level, the magnitudes of /i" {a = 1,2) are of 0{1), by which the cross section 
becomes very large. Furthermore, thanks to the Majorana nature of the t-channel diagram, 
we obtain a much larger cross section in the e~e~ collision than at the e~^e~ collision when 
3> s. The cross section can be as large as 10 pb for m^± = 400 GeV for y/s^-^- = 1 
TeV, m^i = m^2 = 3 TeV and hi = h? = 2.0: see Fig. [HI The backgrounds are expected 

R R 

to be much less than the e^e~ collision. 

We emphasize that a combined study for these processes would be an important 
test for our model, in which neutrino masses are generated at the three loop level 
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by the Z2 symmetry and the TeV-scale RH neutrinos^. In the other radiative seesaw 
models in which the neutrino masses are induced at the one-loop level with RH neu- 
trinos, the corresponding coupling constants to our couplings are necessarily one or 
two orders of magnitude smaller to satisfy the neutrino data, so that the cross section of 
the t-channel RH neutrino mediation processes are small due to the suppression factor (/i" 

Lepton flavor violation: 

Finally, the couplings h'^ cause lepton flavor violation such as /i — > 67, depending on mAr^. 
If such a phenomenon is observed at future experiments [70], we could obtain information 
on m^^. In addition, our model predicts -B(/i — > 67) ^ i?(r — > 67) ^ i?(r — ^ /i7). 

In summary, the possible scenario in this model provides discriminative phenomenological 
characteristics so that it can be tested at future experiments. 



We have discussed the various features of this model neglecting the CP violating phases in 
the Higgs sector, which are crucial for generating baryon number at the EWPT. We comment 
on the case with the CP violating phases. Our model includes the THDM, so that the same 



generation mechanism can be applied in eva. 
at the EWPT unless tan/3 is too large 23 



nation of produced baryon number asymmetry 



24( 1 . For a larger value of tan/5, the constraint 



from the EDM data would be more serious. The mass spectrum in the Higgs sector would 
be changed by including the CP violating phases, but most of the phenomenological features 
discussed above should be conserved with a little modification. 



VII. DISCUSSIONS AND CONCLUSIONS 



We have discussed the model in which neutrino oscillation, DM, and baryon asymmetry 
of the Universe can be simultaneously explained by the TeV-scale physics without intro- 
ducing large fine tuning. Tiny neutrino masses are generated at the three loop level due 
to the exact symmetry, by which stability of the DM candidate is also guaranteed. 
The extra Higgs doublet is required not only for the tiny neutrino masses but also for 

^ Unlike our model, in the model in Ref . |l17| , the coupling constants corresponding to our /i" are small and 
instead those to /i^ are 0(1), so that its Majorana structure is not easy to test at e^e~ collisions. 
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successful electroweak baryogenesis. The phenomenology of the model has been discussed, 
and it has been found that there are several successful scenarios under the constraints from 
the current experimental data. The predictions have been discussed in these scenarios at 
the present and future experiments. It turns out that the model provides discriminative 
predictions especially in Higgs physics and DM physics, so that it is thoroughly testable in 
future experiments. 
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APPENDICES 



A. SM-like limit 

Especially, for the case of the SM like limit (sin(Q; — (3) = —1 |7l|, namely a = /3 — 7r/2, 
(sin a = — cos/5 and cos a = sin/5)), h becomes the SM-like Higgs boson and H decouples 
from the gauge fields (with respect to the three point couplings). In this case, the Yukawa 
interactions are 



/'Quarks 



and 



rrit cot 3 ^- rrit , - 
-nittt H -Htt -htt 

V V 

—mj}b H —Hbb -hbb 

V V 

'rrit -■ 'rrit ^ 
H ztt'y^t H cot pAtt'y^t 

V V 

-zbi'-f^b cot jSAbi^^b 

V V 



+uj'b 
+H-b 



/'Leptons 



\plmt 1 + 75 V2mh 1 — 



75 



V 2 
V2mt cot /5 1 + 75 



t + h.c. 



V 2 

\/2mb cot P 1 



75 



t + h.c. 



(69) 



V2mr _ 1 + 75 + , , 
H u TUJ + h.c. 



V 2 
y/^rrir tan /? _ 1 + 75 



V 2 
mr tan (3 



tH+ + h.c. 



Htt -hTT 



V V 

"rrir _. m-r n a-- 
ZTt'y^T H tan pATt'y^T. 

V V 



(70) 



In the decoupling limit (sin(a — /3) = —1), the masses of scalar fields are expressed as 



ml 



Ml^ = (Ai cos^ f3+X2 sin^ (3 + 2X cos^ (3 sin^ f3)v\ 

= + (-Aicos^/5 + A2sin^/3 + Acos2/3) cos/3sin/3w^ 



2 

Ast;^ 



5 2 



m 



|± = (piCos2/? + p2sin2/5) — , 



/i^ + ((Ti cos^ /? + (T2 sin^ f3)—- 



(71) 
(72) 
(73) 
(74) 
(75) 

(76) 
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B. Neutrino mass matrix 



In order to compute the three loop induced neutrino mass matrix M^j in Eq. fl2^ . we 
start from the calculation of the effective vertex of ij^{pi)-i^{p2)-H^{p3)-H^{p4) comes from 
the loop by Z2-odd particles, iT{pi,p2,P3,P4), where pi, p2, Ps and p^ are incoming momenta 
of Ir^pi), £^(p2), H+{p3), and H^{pa) respectively; 

iTij{pi,p2,P3,PA) 



a=l 



x{~tV2Kv)- ^^ -{-tV2Kv)- ^ (77) 

(ft + pi +P3)^ - [k + pi + ps + p4j^ - mg± 

The Majorana mass matrix of the left handed neutrinos is generated by connecting the 
external lines of Tij{pi,p2,P3,Pi), and with v)^ (also with v)^) by the Yukawa 

coupling, and integrate all the internal momenta taking into account the momentum con- 
servation. There are two Feynman diagrams which give exactly the same contribution, 
corresponding to the way of connecting which couples to v]^: see Fig. [H The mass 
matrix is calculated as 

d^Pi f d^p2 i , ...sMx.., m A -75^ ^(]/i+"^£j 



J i2n)D J {2n)Dpl-ml^^ 2 J pj - 

X iTij {pi , P2 , -Pl , -P2 ) ^ 2^ ^2 0^ 1 i-'^vT^^^f^)— 2 



d^Pl f d^p2 i I SMx„„ A-75^ ^(]/l+"^£j 



(27r)^ y {271)0 p\-m\^^ 'V^. ^^npj ^ ^ ) ^2 _ ^2 
xiTij{pi,p2, -P2, -Pi)— — — - — {-iVe, tanfj) 



2 2 \ o IV till 2 2 

P2 - v 2 y p^- mjj± 



+4/.Vtan^/3 (1^) [(l/r^Dd/r^^) {/ 

d^Pi ]/i 1 1 

(27r)^ - m], p\ - m^± {k + pi)^ - m|± 

d''p2 (-1/2) 1 



(27r)^ _ rn\ja k"^ — rnl 



(78) 



(27r)^ (-^2)^ - "^^j {-P2Y - {k + (-^2))^ - "^|± . 

where we used momentum conservation law ps = —pi and p^ = —p2 {ps = —pi and p^ = —pi) 
for the first (the second) diagram, neglecting the invariant mass of neutrinos. By using 
Passarino-Veltman formalism for one- loop integral functions |33(], we obtain the expression 
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as 



iM, 



+ 



X 



167r2 



1-75 



2 

E 



m 



X [Bi{k'^,mjj±,ml±) - m|, m|±)) 



(79) 



where Bi{k'^, mi, 7712) is the tensor coefficient by Passarino and Veltman 33|. As -C tti^^, 
we neglect in the expression and obtain 



iMij = + 



m 



1 Wl -75\ 4/s:Vtan2/3 ^ 



167r^ 



/d k 
{Bi{k^, m^±, m|±) - m^., m|±)) 

As the asymptotic behavior k^ ^ 00 oi the -Bi function is 



mz 



k"^ — m'jfc, k'^ — 



(80) 



(81) 



so that the integral over d'^k is not divergent. For numerical evaluations, we work in the 
Euclideanized momentum space. 



(82) 



and introduce the cutoff scale A which is a very large number as compared to the scale of 
rrifj or rrig and so on. Then we obtain the expression in Eq. fl24l) . 



C. Line photon flux from DM annihilation 



The differential flux of gamma-ray from a DM annihilation near the center of our galaxy 
is given as 



dn 



{av{-^ ^^))5{E - rrijj) + (erf (- 



dmp{if, (83) 



l.o.s 



where ip is the angle to the galactic center direction, p{l)Js the mass density distribution for 
the DM and we will integrate it along the line of sight / 72| . The first term in the right hand 
side of Eq. ( 183|) denotes the line spectrum comes from the annihilation into two photons, 
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and the second term does to the continuous one. The latter with the differential photon 
spectrum dN/ dE dominantly comes from the decay of pions produced by the fragmentation 
or decay of DM annihilation final state / such as b-quark or r lepton. 

For the integration around the line of sight axis over the solid angle Afi, we introduce a 
dimensionless function \Jm 



in which all informations about the DM halo model are encoded. For a given Af2 ~ 
10~^(10~^), one may find, for instance, J ~ 10^(10^) for NFW density profile [gsI and 
J ~ 10(10) for isothermal 7^. With the integration over angle and substitution of Eq. (18^ . 
we obtain Eq. fl65l) . 



D. High Temperature Expansion 



Let us calculate the integral If by using the high temperature expansion 5J]. For the 
case of bosonic contrbutions, we have 



I Bio) 



dxx^ log[l — e 



27r2 

ji4 



^log[l-e-^1 



where u = y^p^ + m^. The integral is expanded as 
1 r d^p 



f3Jo iM' 



log[l — e 



90" 



^ rr4: ^ 2rr2 , 3rr , ^ 

-I ml H ml + 



log 



24 " ' Utt " ' 647r2 
where 'j^ = 0.5772 is the Euler constant. For fermions, we obtain 

-\/x'^-i-a'^~\ 



+ 



poo 

/ dxx"^ log[l + 
Jo 

27r2 r 1 r d^p 



271^ (771^ 



720 



(27r)3 
1 
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log[l + e-^"] 



rr4 ^ 2rr2 ^ 

I ml — 



647r2 
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m 
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